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ABSTRACT 

Aims. We present the analysis of a local (z = 0.04 — 0.2) sample of 31 galaxy clusters with the aim of measuring 
the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set 
constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. 

Methods. We exploit the large held-of-view and low instrumental background of ROSAT /PSPC to trace the density of 
the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r 200 
and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles 
with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density 
and scatter profiles with the results of numerical simulations. 

Results. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra , we 
observe a steepening of the density profiles beyond ~ r^oo- Comparing our density profiles with simulations, we find 
that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or treating gas 
clumping are in better agreement with the observed gas distribution. We report for the first time the high-confidence 
detection of a systematic difference between cool-core and non-cool core clusters beyond ~ 0.3r2oo, which we explain by 
a different distribution of the gas in the two classes. Beyond ~ rsoo, galaxy clusters deviate significantly from spherical 
symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the 
observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. 
Conclusions. Comparing our results with numerical simulations, we find that non-radiative simulations fail to reproduce 
the gas distribution, even well outside cluster cores. Although their general behavior is in better agreement with the 
observations, simulations including cooling and star formation convert a large amount of gas into stars, which results in 
a low gas fraction with respect to the observations. Consequently, a detailed treatment of gas cooling, star formation, 
AGN feedback, and taking into account gas clumping is required to construct realistic models of cluster outer regions. 
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1. Introduction 


The outskirts of galaxy clusters are the regions where 
the transition between the virialized gas of clusters and 
the accreting matter from large-scale structure occurs and 
where the current activity of structure formation takes 
place. Around the virial radius, the assumption of hydro- 
static equilibrium, which is a necessary assumption for 
the reconstruction of cluster masses from X-ray measure- 
ments, might not be valid any more (e.g., lEvrard et al.1 
1.19961), which could introduce biases in X-ray mass proxies 
dRasia e t al.ll2004: Piffaretti fc Valdarnini [20081: lLau et al.l 
120091: iMeneghetti et al.ll201f)HFabian et al.ll201lli . As a re- 
sult, the characterization of the X-ray emitting gas in the 
outer regions of galaxy clusters is important for mapping 


the gas throughout the entire cluster volume, studying the 
formation processes currently at work in the Universe, and 
performing accurate mass es timates for cosmological pur- 
poses (e.g., lAllen et~ahll2011l) . 

Because of the low surface-brightness of the X-ray 
emitting gas and the extended nature of the sources, 
measuring the state of the int ra-cluster gas around th e 
virial radius is challenging (lEttori fc Molendil 1201 ll) . 
Recently, the Suzaku satellite achieved a breakthrough 
this domain, performing measurements of clus- 
temperatures out to r 20cB 


m 

ter 


. dReiorich et al.l 120091: 

iBautz et al.l [2009: Ka waharada et alT 20101 : Hoshino et, al l 


1 We define ca as the radius within which the mean density 
is A times the critical density 
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201 nt ISimionescu et ahl 
Hum phrey et all 20 1 ill . 
ronn ([George et al l l2009h 


2011 : lAkamatsu et al . 2011 : 
and even in one case beyond 
although the latter detecti on 


is like ly hampered by systematic effects (jEckert et al.1 
l2011aT) . Interestingly, some of the Suzaku results indicate 
very steep temperature profiles and shallow density 
profiles in cluster ou tskirts, at variance with the results 
from XMM-New t on dLeccardi fe; Molendi 2008T) . C handra 


(|Vikhlininet L _al ; 2006' Ettori fc Balestral l2009f) . and 


BOS AT (|Vikhlinin el all 119991: IN euman n| l20Q5|), and with 
the results from numerical s i mulations ( Roncarelli et al.l 
l2006t iTozzi fc Normanl 120011 : iNaeai fc Laul 2011 1. Thus, 
the behavior of the gas in cluster outskirts is still the 
subject of debate. 

Thanks to its large field of view (FOV, ~ 2 deg 2 ) 
and low instrumental background, ROSAT /PSPC is to the 
present day the most sensitive instrument to low surface- 
brightness emission. Its ability to de tect cluster emis s ion at 
larg e radii has been de monstrated bv lYikhlinin et al] Jl999) 
and I Neuman 3 (1200511 (hereafter, V99 and N05). Because 
of the large FOV, it can perform simultaneous local back- 
ground measurements, and therefore it is less affected than 
Suzaku by systematic uncertainties. Its main limitation, 
however, is the restricted band pass and poor spectral res- 
olution, which makes it impossible to measure cluster tem- 
peratures. 


In this paper, we present the analysis of a sample of 
31 galaxy clusters observed with AOSAT/PSPC, with the 
aim of characterizing the cluster emission at large radii and 
comparing the r esults with three different sets of numeri- 
cal simulations (fkoncarelli et al.l [20001 : iNagai fc Laull201lt 
IVazza et al.l l2010l) . The paper is organized as follows. In 
Sect. El we describe our cluster sample and the available 
data. We present our data analysis technique in Sect. [3] and 
report our results in Sect. 0] We compare our results with 
numerical simulations in Sect. [5] and discuss them in Sect. 

Throughout the paper, we assume a ACDM cosmology 
with fl m = 0.3, Ha = 0.7, fib = 0.047 and H 0 = 70 km s -1 
Mpc -1 . 


2. The sample 

We select objects in the redshift range 0.04 — 0.2, such that 
r 2 oo is easily contained within the FOV of the instrument 
and is large enough to allow for a sufficient sampling of 
the density profile. We restrict ourselves to observations 
with sufficient statistics to constrain the emission around 
the virial radius. Our final sample comprises 31 clusters 
in the temperature range 2.5-9 keV, with the addition of 
A2163 ( kT ~ 18 keV). Among our sample, we classify 
14 clusters as cool - core ( CC) following the classification 
of ICavagnolo et aD (120091 1 (i.e. they exhibit a central en- 
tropy I\ 0 < 30 keV cm 2 ), and 17 as non-cool core (NCC, 
A' 0 > 30 keV cm 2 ). We recall that CC clusters exhibit a 
relaxed morphology, a high central density and a temper- 
ature decrement in the central regions, while NCCs trace 
dynamically-disturbed clusters with irregular morphologies 
and flat temperature and density profiles in their cores (e.g., 
ISanderson et alJfeOOOl : iHudson et al.ll2010ll . 

Our sample of clusters, together with the log of the 
available data and some important quantities, is shown in 
Table [3] In Fig. [T] we plot the distribution of temperature 


(left panel) and central entropy (right panel) for our sam- 
ple. It should be noted that the sample was selected based 
on the quality of the existing observations, and might be 
subject to selection effects. However, for the purpose of this 
work we do not require that the sample be representative or 
complete, since we are interested in the characterization of 
cluster outskirts, which exhibit a high level of self-similarity. 


3. Data analysis 

3.1. Data reduction 

We use the ROSAT Extended Source Analysis Software 
(jSnowden et al.lfl994h for data reduction. We filter out time 
periods when the master veto count rate exceeds 220 cts/sec 
(using valid_times), and extract light curves for the whole 
observation using rate_pspc. We use the ao executable to 
model the atmospheric column density for the scattering of 
solar X-rays, and fit the light curves in each energy band to 
get the relative contributions of the scattered solar X-rays 
(SSX) and of the long-term enhancements (LTE) , using the 
rate_fit executable. 

We then extract event images in each energy band 
and the corresponding effective exposure maps, taking 
into account vignetting effects. We compute the contri- 
bution of the various background components, the LTE 
(lte_pspc), the particle background (cast_part), and the 
SSX (cast_ssx), and combine them to get a total non- 
cosmic background map. 

3.2. Surface-brightness profiles 

The point-spread function (PSF) of ROSAT/ PSPC is 
strongly angle- dependent, and ranges from ~ 15 arcsec on- 
axis to 2 arcmin in the outer parts of the FOV. Thus, the 
sensitivity of the instrument to point sources is higher on- 
axis, and a larger fraction of the cosmic X-ray background 
(CXB) is resolved. Consequently, when detecting sources 
in the image it is important to use a constant flux thresh- 
old, such that the same fraction of the CXB is resolved all 
over the FOV and the value measured in the source-free 
regions can be used to subtract the background. We detect 
point sources using the program detect with a minimum 
count rate of 0.003 cts/sec in the R3-7 band (~3x 10“ 14 
ergs cm' 2 s _1 in the 0. 5-2.0 keV band) to resolve the same 
fraction of the CXB over the FOV, and mask the corre- 
sponding areas. To compute surface-brightness profiles, we 
extract count profiles from the event images in the R3-7 
band (0.42-2.01 keV) with 30 arcsec bins centered on the 
surface-brightness peak, out to the radius of 50 arcmin. 
We divide each pixel by its corresponding exposure to ac- 
count for_the vignetting effects, following the procedure of 
lEckert, et al.1 (1201 lliPl . We perform the same operation for 
the background map and subtract the non-cosmic back- 
ground profile in each bin. 

We tested this procedure on 4 different blank fields to 
estimate the accuracy in our determination of the CXB. 
We extracted the surface-brightness profile for the 4 obser- 
vations from the center of the FOV, grouped the bins to 
ensure a minimum of 100 counts per bin, and fitted the 
resulting profiles with a constant (see Fig. El). While the 
agreement is qualitatively good, significant deviations to 

2 http:/ /www. iasf-milano.inaf.it/~eckert/newsite/Proffit.html 
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Fig. 1. Distribution of temperature (left) and central entropy (right) of the members of our sample (see Table [3J. In the 
left panel, A2163 ( kT ~ 18 keV) is located outside of the range. 



Fig. 2. Surface-brightness profiles for 4 blank-field PSPC 
observations from the center of the FOV, fitted with a con- 
stant. The dashed line shows the vignetting correction curve 
for comparison, in arbitrary units; the bump at ~ 22 ar- 
cmin is caused by the support structure. The bottom panel 
shows the ratio between data and model. 


the model are found, leading to an excess scatter of ~ 6%, 
which we use as an estimate of the systematic uncertainties 
in the measurement of the CXB. This value encompasses 
both the cosmic variance and the true systematic uncer- 
tainties, e.g., in the vignetting correction or the determina- 
tion of the particle background. The higher level of scatter 
in the central regions is explained by the small area of the 
corresponding annuli, which implies a large cosmic variance 
likely due to discrete sources with fluxes just under our ex- 
clusion threshold. Since, in most cases, the value of r 2 oo is 
larger than 15 arcmin, our systematic error of 6% is a con- 
servative estimate of the level of systematic uncertainties 
at the virial radius. 

For each cluster, we then use te mperature pro- 
files from the literature ( XMM-Newton , ISnowden et al.l 


J2008 ) ; Chand ra , Cavagnolo et al.1 (l2009h ; BeppoSAX , 
iDe Grandi fe Molendil ( 2002f) j to estimate the virial tem- 
perature of the cluster. We approximate T v i r as the mean 
temperature in the 200-500 kpc region, i.e. excluding the 
cool core and the temper ature decline in the outskirts 
(lLeccardi fc Molendil 12008 ) . Using this estimate of T vir , 
we comp ute the val ue of r 2 oo from the scaling relations of 
lArnaud et al.l (1200511 . We then use the source- free region of 
the observation (r > 1.3r 2 oo) to fit the surface-brightness 
profile with a constant and get the cosmic background level 
for the observation, with the exception of the Triangulum 
Australis cluster, for which we use the range r > l.lr 2 oo 
because of the large value of r 2 oo (~ 37 arcmin). 

After having estimated the sky background for our ob- 
servation, we extract again the surface-brightness profile in 
the radial range 0 — 1.3r2oo with logarithmic bin size. The 
best-fit value for the CXB is subtracted from the profile and 
its error is added in quadrature to each bin. The systematic 
error of 6% on the CXB is also added in quadrature to ac- 
count for the cosmic variance and systematic uncertainties. 
For comparison, we note that in most cases the statistical 
uncertainties in the profiles are of the order of 10% of the 
CXB value around r 2 oo- 


3.3. Density profiles 

To compute the density profiles, we first rebin our 
background-subtracted surface-brightness profiles to ensure 
a minimum of 200 counts per bin and a detection signifi- 
cance of at least 3 <t, to reach suffi cient s tatisti cs in each bin. 
We use the procedure of iKriss et al.l (1983) to deproject 
the observed profiles, and the PSPC response to convert 
the observed count rates into emission measure, through 
the no rmalization of the MEKAL model (see lEckert et al.l 
l2011ai for details), 

10 -14 f 

Narm = toK,(i + s)p J n ‘ nHdV ’ (1) 

which is proportional to the emission measure. To this aim, 
we assume that the spectrum of our sources is described 
by an absorbed MEK AL model with Njj fixed to the 21cm 
value dKalberla et al.ll2005lf . abundance fixed to 0.3Z Q , and 
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the temperature profiles adopted from the literature (see 
Table [3]), and fold the model with the PSPC response. The 
conversion from PSPC R3-7 count rate to emission measure 
is then inferred. Beyond the limit of the temperature pro- 
files, the temperature of the outermost annulus is used. We 
note that the conversion from PSPC count rate to emission 
measure is largely insensitive to the temperature: between 
2 and 8 keV the conversion factor changes at most by 4%. 
Once converted into the MEKAL normalization, we infer 
the density profiles, assuming spherical symmetry and con- 
stant density into each shell. 

The error bars on the density profiles were estimated us- 
ing a Monte Carlo approach. In every case, we generated 10 4 
realizations of the surface-brightness profile using Poisson 
statistics, and performed the geometrical deprojection fol- 
lowing the method described above. The la error bars were 
then estimated by computing the root-mean square devia- 
tion (RMS) of our 10 4 realizations of the density profile in 
each density bin. 


3.4. Azimuthal scatter profiles 

For the purpose of this work, we are also interested in the 
deviations of the X-ray emission from spherical symmetry. 
To this aim, we divide our images into N azimuthal sec- 
tors with constant opening angle, and compute the surface- 
brightness profiles in each sector individually. We then 
compute the scatter of the various sectors with respect 
to the mean profile, following the definition introduced by 
IVazza et al.l (1201 lbl f. 


, 2 1 y (SBj - {SB)) 2 

N ^ (SB) 2 

1=1 ' 


( 2 ) 


where (SB) is the mean surface-brightness and SBi,i = 
1..N denotes the surface-brightness computed in the var- 
ious sectors. Since the statistical fluctuations of the data 
also introduce a certain level of scatter, it must be noted 
that the quantity computed through expression [2] gives the 
sum of the statistical and intrinsic scatter, 

E 2 = E 2 „ t + E 2 tat . (3) 

The statistical scatter E stat is given by the mean of the 
individual relative errors, 


1 iV 2 

y2 _ J_ Y' 

stat N ^ (SB) 2 ’ 


( 4 ) 


and must be subtracted from Eq. [2] to estimate the level 
of intrinsic scatter. The validity of the aforementioned for- 
mula for the statistical scatter was verified through a set of 
simulations of a source with no intrinsic scatter. The un- 
certainties in the scatter are then estimated through Monte 
Carlo simulations. 

In our analysis, we group the bins of the total surface- 
brightness profiles to reach a minimum of 8 <7 per bin, and 
then divide our images into 12 sectors with an opening of 
30°. The result of this analysis is a profile describing the 
intrinsic azimuthal scatter of the X-ray surface brightness, 
in percent. 

It must be noted that the method presented here is sen- 
sitive to all kinds of deviations from spherical symmetry, 


whether it is induced by the asymmetry of the large-scale 
structure (e.g., filaments), by gas clumping or by ellipticity. 
The cause of the observed asymmetry cannot be determined 
from the azimuthal scatter alone. 


4. Results 

4.1. Emission measure and density profiles 

In Fig.[3]we show the scaled emission measure profiles (left, 
following Eq. [[]) and the deprojected density profiles (right) 
for the 31 clusters of our sample. A self-similar scaling 
was a pplied to the emission-measure profiles (lArnaud et al.l 
l2002f) . i.e. each profile was rescaled by the quantity 




1/2 


( 5 ) 


The density profiles were rescaled by E 2 (z ) = fl m (l + 
+ Ha follo wing t heir expected evolution with redshift 


(ICroston et alJ 12008ft. As a l ready noted by s everal authors 
(e.g., ViMrlini^etMj 1999|; iNeumannl 120051 ICroston et, al.l 
12008c Leccardi et al.l l2010f) . the profiles show a remarkable 
level of self-similarity outside of the core (r > 0.2r 2O o). 
On the other hand, the large scatter observed in the cen- 
tral regions reflects the distinction of the cluster popu- 
lation into CCs, showing a prominent surface-brightness 
peak, and NCCs, which exhibit a flat surface brightness pro- 
file in their cores, as expecte d from the standard /3- model 
(jCavaliere fc Fusco-Femianolll976l f. 


SB(r) = SB 0 




(6) 


In the radial range 0.2 — 0.7r 2O o, the scatter of the den- 
sity profil es is 10-20%, in exce llent agreement with the 
Chan dra (IVikhli nin et all l2006f) and XMM-Newton re- 
sults (ICroston et al.ll2008h . However, ICroston et al.l (|2008fl 
needed to rescale the profiles by T -1 / 2 to account for the 
lower gas fraction in low-mass objects. In our case, perform- 
ing such a scaling does not further reduce the scatter of the 
profiles. This is probably explained by the relatively nar- 
row temperature range spanned in our sample (all but one 
objects have a temperature higher than 3 keV), such that 
the clusters in our sample should show little dependence on 
gas fraction. 


4.2. Stacked emission-measure profiles 

To compute the mean profile of our sample, we interpolated 
each profile following a pre-defined binning, and performed 
a weighted mean to compute stacked profiles. The errors 
on the interpolated points were propagated to the stacked 
profiles. We also divided our sample into the two classes 
(CC and NCC) to look for differences between them. 

In Fig. [I] we show the stacked emission-measure (EM) 
profile for the entire sample (black) compared to the pro- 
files stacked for the two populations separately (see also 
Appendix [B]) . Interestingly, we note a clear distinction be- 
tween the two classes in cluster outskirts (see the bottom 
panel of the Figure). Namely, beyond ~ 0.3r 2 oo NCC pro- 
files systematically exceed CCs. We stress that this effect is 


D. Eckert et al.: The gas distribution in galaxy cluster outer regions 


5 




Fig. 3. Scaled emission measure (left, in units of cm 6 Mpc) and density profiles (right) for the 31 clusters of our sample 
(see Table©. 



Fig. 4. Stacked emission measure profile (in units of cm 6 
Mpc) for the entire sample (black), and the two populations 
individually (CC, red; NCC, blue). See also Appendix [B] 
The bottom panel shows the ratio between the CC and 
NCC populations. 


really a difference between the two classes, i.e. it is not in- 
troduced by a biased distribution of another quantity (such 
as temperature or redshift). Indeed, grouping the profiles 
according to the temperature or the redshift did not show 
any particular behavior, which indicates that we are really 
finding an intrinsic difference between the CC and NCC 
classes. This result could follow from a different distribu- 
tion of the gas in the two populations or from a higher 
clumping factor in disturbed objects (see Sect.©. 

Alternatively, the observed difference could be explained 
by an inaccurate determination of r? on f or NCC c lusters. 
Indeed, the scaling relations of lArnaud et al. (i2005ll were 
computed under the assumption of hydrostatic equilib- 
rium, which is better fulfilled in CC clusters. This ex- 
planation is, however, unlikely. Indeed, to recover self- 
similarity, our value of r 2 oo should have been system- 
atically under- estimated by ~ 10% for NCCs, i.e. since 


1 /2 

r 2 oo °c T V ir the virial temperature of the NCC clusters 
should have been under-estimated by more than 20%. From 
mock C handra o b seryations o f a sample of simulated galaxy 
clusters. iNagai et~aD (2007b) determined that the spectro- 
scopic temperatures of unrelaxed clusters differs from that 
of relaxed clusters by ~ 5%, which is insufficient to explain 
the observed difference. It is therefore unlikely that such a 
large error on the virial temperature would be made. 

We fitted the mean scaled emission-measure profiles 
from Fig. [I] with the standard /3-model (Eq. ©, adding 
a second /3 component in the case of the CC clusters to 
take the cool core into account. The (double) /3 model 
gives a good representation of the data in the radial range 
0 — 0.7r2oo (~ rsoo), but significantly exceeds the observed 
profiles above thi s radius, in agreement w ith the results 
of V99, N05 and lEttori fc Balestra] (1200911 . For CC clus- 
ters, the best-fit model gives /3 = 0.717 ± 0.005, while for 
NCC clusters we find (3 = 0.677 ± 0.002. Fitting the radial 
profiles in the range 0.65-1.3r2oo, we observe a significant 
steepening, with a slope /3 = 0.963 ± 0.054 for CCs and 
/3 = 0.822 ± 0.029 for NCCs. As explained above, the slope 
of the NCC profile is flatter than that of the CC profile 
beyond r$ oo- In more detail, the fits of the profiles in var- 
ious radial ranges are reported in Table © to quantify the 
steepening. 

Given the limited number of objects in our sample, we 
have to verify that this result is not a chance realization. 
To this aim, we fitted all the emission-measure profiles at 
r > 0.3r 2 oo with a /3 profile, fixing the value of /? to 0.7 and 
r c to 0 . 12 r 2 och and extracted the best-fit normalization for 
all profiles. We then sorted the normalization values into 
the CC and NCC classes, and performed a Kolmogorov- 
Smirnov test to determine the probability that they origi- 
nate from the same parent distribution. Using this proce- 
dure, we found that the chance probability for this result is 
very small, P ~ 6 x 10 -7 . Therefore, we can conclude with 
good confidence that we are indeed finding an intrinsic dif- 
ference between the two classes. 

4.3. Stacked density profiles 

We stacked the density profiles shown in the right panel 
of Fig. © following the same method as the EM profiles. 
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Fig. 5. Average proton density profile for the entire sample. 
The dashed lines indicate the positive and negative scatter 
of the profiles around the mean value. 


From the different profiles, we also computed the scat- 
ter of the profiles around the mean value, following a 
method similar to the one presented in Sect. 13.41 for the 
azimuthal scatter. The statistical scatter was subtracted 
from the total scatter using the same technique. In Fig. 
[5] we show the average density profile of our clusters to- 
gether with the scatter of the individual profiles around 
the mean value (see also Table [TJ. At r 20 o, the mean den- 
sity is 7T2oo = (3.8 ± 0.4) x 10 ~ 5 E 2 (z) cm -3 , with 25% 
scatter. For comparison, it is interesting to note that the 
den sity of PKS_0745-191 claimed in the Suzaku analysis of 
I George et al.l (2009) at r 20 o deviates from our mean value 
by more than Jht, which c asts ev en more doubts on this 
measurement (lEckert et al. 11201 fal l. 

As for the EM, we also extracted mean density profiles 
individually for the two classes of clusters in our sample. 
The same behavior is observed at large radii, i.e. the den- 
sity of NCC clusters is systematically higher (by ~15%) 
compared to CCs above r ~ 0.3r 2 oo- A global steepening of 
the density profiles is also observed beyond ~ rsoo- 

Our density profiles are in good agreement with the re- 
sults of V99. However, while V99 estimated the density 
from /3-model fitting, we performed a geometrical depro- 
jection of the data using temperature profiles to infer the 
mean density profile. This method has the advantage of 
being model-independent. 


4.4. Gas mass 

We computed the gas mass from our deprojected density 
profiles and stacked them in the same way as described 
above. In the self-similar model, the gas mass is expected 
to fo llow the relation M oc T 3 / 2 (e.g., iBrvan fc Norman! 
I1998D . However, observational works indicate that the ac- 
tual M gas — T relat ion is steep er than the expected se lf- 
similar scaling ^N euman n fc Arnaucfl 12001b lArnaud et al.1 
l2007t ICroston et a,l.ll2008l) because of the lower gas fraction 
in groups and poor clusters. For this work, we use the rela- 


Table 1. Mean emission- measure and density profiles com- 
puted from our sample. Column description. 1 and 2: Inner 
and outer bin radius in units of r 2 oo; 3: Emission measure 
rescaled by A ssc in units of cm -6 Mpc; 4: Average pro- 
ton density in units of 10~ 3 cm -3 ; 5: Scatter of the various 
profiles relative to the mean value in percent. 


Rin 

Rout 

ScEM 

n H E(z ) * 

(J 

0 

0.03 

(1.57 ±0.01) x lO’ 8 

12.17 ±0.04 

56 

0.03 

0.06 

(1.086 ± 0.004) x 10” 5 

6.05 ±0.02 

47 

0.06 

0.09 

(6.52 ±0.03) x 10' 6 

3.44 ±0.01 

36 

0.09 

0.12 

(4.46 ± 0.02) x 10' 6 

2.21 ±0.01 

25 

0.12 

0.15 

(3.01 ±0.01) x 10' 6 

1.58 ±0.007 

21 

0.15 

0.18 

(2.09 ±0.01) x 10~ 6 

1.175 ±0.006 

16 

0.18 

0.21 

(1.42 ±0.007) x 10 -6 

0.919 ±0.006 

17 

0.21 

0.24 

(1.05 ±0.006) x 10” 6 

0.728 ± 0.006 

13 

0.24 

0.27 

(8.16 ±0.05) x 10~ 7 

0.598 ± 0.005 

13 

0.27 

0.3 

(6.05 ± 0.04) x 10~ 7 

0.501 ± 0.005 

15 

0.3 

0.33 

(4.64 ± 0.04) x 10~ 7 

0.422 ± 0.004 

13 

0.33 

0.37 

(3.54 ±0.03) x 10~ 7 

0.350 ± 0.004 

16 

0.37 

0.42 

(2.56 ±0.02) x 10~ 7 

0.284 ± 0.004 

14 

0.42 

0.47 

(1.77 ±0.02) x 10~ 7 

0.225 ± 0.003 

14 

0.47 

0.52 

(1.27 ±0.01) x 10~ 7 

0.188 ±0.003 

20 

0.52 

0.59 

(8.58 ±0.12) x 10~ 8 

0.145 ± 0.003 

19 

0.59 

0.66 

(5.74 ±0.10) x 10 -8 

0.120 ±0.003 

18 

0.66 

0.74 

(4.13 ±0.09) x 10 -8 

0.099 ± 0.002 

13 

0.74 

0.83 

(2.66 ± 0.08) x 10' 8 

0.073 ± 0.002 

26 

0.83 

0.93 

(1.59 ±0.07) x 10' 8 

0.059 ± 0.002 

27 

0.93 

1.05 

(8.08 ± 0.59) x 10~ 9 

0.039 ± 0.002 

26 

1.05 

1.17 

(4.75 ± 0.53) x 10' 9 

0.028 ± 0.002 

36 



Fig. 6. Enclosed gas mass profiles for CC (red) and NCC 
systems (blue) . The data were rescal ed b y E(g)fcT~ 1 986 as 
observed in the REXCESS sample (iCroston et al.l 1200 8). 

tion determined from the REXCESS sample (iCroston et al.l 
120081 ) to rescale our gas mass profiles, 

/ urn \ i-986 

«-«£(,)-■(_) . (7, 

As above, we divided the sample into CC and NCC 
classes, and stacked the two classes individually. In Fig. [6] 
we show the mean gas mass profiles for CC (red) and NCC 
clusters (blue). As expected, CCs have a higher gas mass in 
their inner regions, since their central densities are higher. 
More interestingly, we see that the two profiles converge in 
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cluster outskirts, and exhibit a gas mass around the virial 
radius that is consistent within the error bars. At r 2 oo, the 
universal gas mass is 


M gas , 2oo = (2.41±0.05) x \Q 1A E{z)~ 1 


kT 


10 keV 


1.968 


M 0 , (8) 


with a scatter of 17% around the mean value. This result 
follows from the higher density measured in average beyond 
~ 0.3r 2 oo in NCC clusters and the steeper slope of CC pro- 
files in the outskirts (see Sect. 14.21) . The lower density of 
CC clusters in the outer regions compensates for the well- 
known excess observed in the cores, such that the total gas 
mass contained within the dark-matter halo follows a uni- 
versal relation. We also estimated the average gas fraction 
by computing th e expecte d value of M 2 oo using the scaling 
relations of lArnaud et al.l (2005). For our sample, we find 
a mean gas fraction within r 2 oo of 


/ hj' \ 0 478 

= (0.15 ± 0.01) (jj^j 


( 9 ) 


in good agr e ement with previous work s (e.g., 
IVikhlinin et al.l 120061 : iMcCarthv et al.l l2007f) , which 
for the highest mass ob jects corresp o nds t o ~ 89% of the 
cosmic baryon fraction llarosik et alJ[20T ll. 


4.5. Azimuthal scatter 

Following the method described in Sect. 13.41 we computed 
the azimuthal scatter of the surface-brightness profiles for 
all the clusters in our sample, and rescaled the scatter pro- 
files by our estimated value of r 2 oo- We then stacked the 
profiles using the same procedure as described above and 
computed the mean azimuthal scatter. We recall that since 
the surface brightness depends on n \ the variations in den- 
sity are less important than the ones computed here. 

In Fig. □ we plot the average scatter profile (black), 
compared to the mean value for CC (red) and NCC clus- 
ters (blue). The increase in the innermost bin is an artifact 
introduced by the small number of pixels in the center of 
the images, and therefore it should be neglected. At small 
radii (r < 0.5r 2 oo) we find a clear difference between CC 
and NCC clusters, that is easily explained by the more dis- 
turbed morphology of the latter. In this radial range, CC 
profiles exhibit a scatter of 20-30%, which corresponds to 
density variations of the order of 10% , in good agreement 
with the value predicted by lYazza etah (201Tb) from nu- 
merical simulations. Conversely, beyond r ~ r 500 the pro- 
files for CC and NCC clusters are similar, and indicate a 
large scatter value (60-80%). 

We investigated whether any systematic effect could af- 
fect our result in cluster outskirts, where the background 
is dominating with respect to the source. Indeed, in such 
conditions, the total scatter is dominated by the statistical 
scatter. In case the mean level of systematic uncertainties in 
the CXB reconstruction exceeds our adopted value of 6%, 
Eq. [3] immediately implies that the intrinsic scatter would 
be over-estimated. The presence of both intrinsic and sta- 
tistical scatter could also introduce some covariance term, 
which is not taken into account in Eq. [3] To test this hy- 
pothesis, we ran a set of simulations including source and 
background, where we introduced a given level of intrin- 
sic scatter for the source and a systematic error in addi- 
tion to the Poisson statistics for the background. We then 



Fig. 7. Stacked azimuthal scatter profile for the entire clus- 
ter sample (black). The red and blue data represent the 
mean profile extracted from CC and NCC clusters, respec- 
tively. 


computed the intrinsic level of scatter following Eq. [3] Our 
simulations indicate that even when increasing the level of 
systematic uncertainties to 12% of the CXB value, a sig- 
nificant bias in the measured scatter only appears when 
the source-to-background ratio is of the same order as the 
systematic uncertainties. Since, by construction, we never 
detect a signal when the source is less than ~15% of the 
CXB value, our results are unaffected by these effects, and 
we can conclude with good confidence that the high level 
of scatter measured beyond ~ rsoo is an intrinsic property 
of our cluster sample. 

V99 also investigated the deviations from spherical sym- 
metry by measuring the value of /? in 6 sectors in the radial 
range r > 0.3rigo, and concluded that the assumption of 
spherical symmetry is relatively well satisfied in cluster out- 
skirts, at variance with our results (see Fig. 0 . However, 
when fitting a /3-model the fit is mostly driven by the shape 
of the profile in the innermost region, where the statistics is 
higher. Conversely, our method is model-independent, and 
directly stacks the data at similar radii. For relaxed objects, 
our data also indicate little deviation from spherical sym- 
metry at r < r.500, and a significant scatter is only observed 
beyond rsoo, so it is probable that these deviations would 
not be reflected in the /3-model fit. For instance, the case 
of A2029 is striking. While, in agreement with V99, we find 
little azimuthal variations of f3 outer, we observe a high level 
of scatter in this object beyond r.soo, which is explained by 
the presence of a possible filament connecting A202 9 to its 
neighbor A2033 in the North (see iGastaldello et all (|2010H 
and Appendix lAl) . Moreover, V99 deliberately excluded a 
number of systems with obviously disturbed morphologies, 
such as A3558 and A3266, which we included in our sam- 
ple. Therefore, our results are not in contradiction with the 
ones of V99. 


5. Comparison with numerical simulations 

In this section, we compare our observational re- 
sults with three differ e nt sets of numerica l simulations 
dRoncarelli et al.1 120061 : iNaeai et al.1 l2007bt lYazza et al.l 
l2010h . We analyze the results of a composite set of cosmo- 
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logical runs, obtained by the different authors with slightly 
different cosmological and numerical setups. In addition, 
the preliminary data reduction was made on each data-set 
following independent post-processing techniques, aimed at 
assessing the role of gas clumping on the comparison be- 
tween simulated mock and real X-ray observations. Our 
aim in this project is to test the most general and converg- 
ing findings of such different runs, against our observations 
with ROSAT/ PSPC. 


5.1. Simulations 

5.1.1. ENZO 

We use a sample of 20 simulated clusters from the high 
resolution and non-r adiat i ve res imulations of massive 
systems presented in [Vazza et al.l (|2010h . In this set of 
simul ations, Adaptive Me sh Refinement in the ENZO 1.5 
code dNorman et al.1 f2007h have been tailored to achieve 
high resolution in the innermost regions of clusters (follow- 
ing the raise of gas and DM overdensity), and also in the 
outermost cluster regions, following the sharp fluctuations 
of the velocity field, associated with shocks and turbulent 
motions in the icm For a detailed presentation of the 
statistical properties of the thermal gas (and of turbulent 
motions) in t hese s imulat ed systems we refer the reader to 
IVa.zza et al.l <1201 ft 1201 lah . 


5.1.2. ART 

We analyze a sample of 10 simulated cl usters with 
Tx > 2.5 keV from the sample presented in iNagai et al.1 
(l2007atf blh These simulations are performed using the 
Adapt ive Refinement Tree (ART) N-body+ gas-dynamics 
code dKravtsovl Il999t iKravtsov et, al.l l2002h , which is an 
Eulerian code that uses adaptive refinement to achieve 
high-spatial resolution (a few kpc) in self-consistent cosmo- 
logical simulations. To assess the impact of cluster physics 
on the ICM properties, we compare two sets of clusters sim- 
ulated with the same initial conditions but with different 
prescription of gas physics. In the first set, we performed hy- 
drodynamical cluster simulations without gas cooling and 
star formation. We refer this set of clusters as non-radiative 
(NR) clusters. In the second set, we turn on the physics of 
galaxy formation, such as metallicity-dependent radiative 
cooling, star formation, supernova feedback and a uniform 
UV background. We refer this set of clusters as cooling+star 
formation (CSF) clusters. For detailed descriptions of the 
gas physics a nd mo ck X-ray ima ges we refer the read er to 
INagai et al.1 (l2007al lbh. Following INagai fc Lad (1201 lh . we 
compute the dumping-corrected gas density profiles of X- 
ray emitting gas with T > 10 6 K for comparisons with 
X-ray observations. Indeed, the formation of dense clumps 
increases the emissivity of the gas, which leads to an overes- 
timation of the measured gas density when the assumption 
of constant density in each shell is made. 


3 A public archive of the final output of these simulations is 
available at http://data.cineca.it 


5.1.3. GADGET 

This set includes 4 massive h alos sim ulated with the 
GADGET-2 Tree-SPH code (ISpringell l2005f). with M 2 on > 
IOWMq (for a detailed description see iRoncarelli et al.l 
(i2006i f and references therein). Each object was simu- 
lated following two different physical prescr i ptions: a non- 
radiative run (referred to as ovisc in IRoncarelli et al.ll2006l f 
and a run including cooling, star formation and supernovae 
feedback (CSF). 

In order to eliminate the dense clumps that dominate 
the density and surface brightness in the outskirts, when 
computing the profiles for every radial bin we excise the 
1 per cent of the volume that corresponds to the densest 
SPH particles. This empirical method mimics the proce- 
dure of masking bright isolated regions from the analysis of 
observed clusters. 


5.2. Comparison of gas density profiles 


We compared the simulations with our observed mean 
ROSAT density profile (see Fig. [5] and Table [T}. We present 
the detailed comparison in Fig. [8} with the non-radiative 
(NR) simulations (left panel) and with the CSF simulations 
(right). From the figures, we find a relatively good agree- 
ment between all the different sets of simulations, espe- 
cially beyond ~ 0.7r2oo- The non-radiative GADGET run has 
a lower normalization than the corresponding grid codes, 
because in GADGET the fraction of baryons virializing into 
clusters is smaller than the cosmic value (~ 78% of the 
cosmic baryon fraction), while grid codes predict a baryon 
fraction in clusters very close to the cosmic value. In gen- 
eral, we see that the predicted density profiles are too steep 
compared to the data. We note that NR runs (ENZO, red; 
ART, dotted cyan; GADGET, dashed green) predict steeper 
profiles than the runs including cooling, star formation and 
feedback effects (ART, magenta; GADGET, dashed blue). CSF 
profiles also have lower normalizations, since radiative cool- 
ing transforms a fraction of the gas into stars. The profile 
including the effects of clumping (dotted magenta) shows 
the best agreement with the data. 

To quantify this effect, we fitted the various pro- 
files in three different radial ranges (0.2 — 0.4r2och 0.4 — 
0.65r2oo, and 0.65 — 1.2r2oo)- In the inner regions, the 
effects of additional physics are expected to be impor- 
tant, thus highlighting the differences between NR and 
CSF runs. The radial range 0.4 — 0.65r2oo(~ 0.6 — R'500) 
is a good range for the comparison with the data, since 
the effects of radiative cooling should be small, and 
data from several different satellites are available for 
cross-check. On the observational side, the density pro- 
files in this radial range are well-fitted by the /3-model 
(see Eq. [6]), and several independent works con ve rge to 
the canonical value of B ~ 0.7 (e.g., iMohr et al.l 119991: 
Ettori&JAbian 19991: V~khlinin et al.ll 199f)t Crost on et all 


20081 : Ettori fc Balestral 20091 : I Eckert et al.l 120 1 111) . As a 
benchmark, we computed the values of /3 for our average 
density profile and the various sets of sim ul atio ns, fixing the 
core radius to 0.12r2oo fe.g.. IMohr et al.lfl999if . The results 
of this analysis are shown in Table [2] The fits to the ob- 
servational data were performed on the emission-measure 
profiles (see Sect. 14.211 to take advantage of the larger num- 
ber of bins and minimize the uncertainties linked to the 
deprojection procedure. 
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Fig. 8. Comparison between the mean ROSAT density profile for our sample and the different sets of numerical simula- 
tions. The shaded area indicates the data and ler scatter as shown in Fig. [5] The bottom panels show the ratio between 
data and simulations as a function of radius. Left: Comparison with non-radiative simulations. The dotted red cu rve 
represents the ENZO profile (IVazza et alJ l201(jh . the solid green curve shows the ART simulations (INagai et al.ll2007bh . 
and the dashed blue curve is the GADGET profile (jRoncarelli et al.l [2006b Right: Same with CSF simulations. The dashed 
blue line shows the GADGET sim ulations, w hile the green curves show the ART profiles, for the total density (solid) and 
corrected for clumping (dotted. iNagai fc Laull201llh 


Table 2. Values of the /? parameter 

(ICavaliere fe Fusco-Femianol 11976T ) in several radial 
ranges for the average ROSAT profiles and the various 
sets of simulations. The core radius was fixed to 0.12r2oo in 
all cases. The subscript cl indicates the profiles corrected 
for the e ffect of clu mping using the method described in 
iNagai fc haul (1201 ill . 


Data set 

do. 2 - 0 . 4 

do. 4-0. 65 

do. 65-1. 2 

Data, total 

0.663 ± 0.002 

0.720 ± 0.009 

0.886 ± 0.024 

Data, CC 

0.702 ± 0.004 

0.712 ±0.016 

0.979 ±0.053 

Data, NCC 

0.638 ± 0.003 

0.731 ±0.011 

0.855 ±0.029 

ENZO 

0.744 

0.945 

0.952 

ART, NR 

0.801 

0.956 

0.983 

ART, CSF 

0.808 

0.842 

1.005 

ART, NR, cl 

0.701 

0.824 

0.854 

ART, CSF, cl 

0.803 

0.718 

0.902 

GADGET, NR 

0.856 

0.857 

0.971 

GADGET, CSF 

0.756 

0.864 

0.944 


These numbers confirm the visual impression that the 
simulated gas density profiles are steeper than the observed 
ones. In the 0.4 — 0.65r 2 oo range, while all our datasets 
converge to a /3 value very close to the canonical value, 
all the simulations lead to significantly steeper gas profiles, 
with /? values larger than 0.85, with the exception of the 
ART profile including CSF and clumping. Therefore, we can 
see that at this level of precision the effects of additional 
physics cannot be neglected, even in regions well outside of 
the cluster core. 

The results presented in Table [2] also highlight the dif- 
ferences between NR and CSF runs. Inside r 5 oo, the simu- 
lations including additional physics lead to flatter density 


profiles compared to the NR runs. In this case, gas cool- 
ing is converting a fraction of the X-ray emitting gas into 
stars. Since the cooling efficiency decreases with radius, 
more gas disappears from the X-ray range in the central 
regions, which results in flatter density profiles and lower 
normalizations. We note, however, that this effect is prob- 
ably overestimated in the CSF simulations. Indeed, it is 
well-known that these simulations predict a stell ar fracti on 
which is well above the observed value ('e.g. JKravtsov et al.l 
120051 : iBorgani fc Kravtso\il200f)!l . This effect is particularly 
strong in the ART CSF simulation, for which nearly one 
third of the gas is converted into stars. Beyond r 50 o, there 
is little difference between NR and CSF runs, i.e. the effects 
of additional physics are not i mp orta nt. At large radii, the 
effect of gas clumping (INagai fe Laull201lli dominates and 
flattens the observed profiles. As we can see in Table [2] 
and in the right panel of Fig. [5J the ART profile including 
both additional physics and a post-processing treatment of 
clumping reproduces better the behavior of the data, even 
though it is still slightly too steep. 

5.3. Azimuthal scatter 

A study of the azimuthal scatter in the radial profiles of 
density, temperature, entropy and X-ray bright ness o f sim- 
ulated ENZO clusters has been presented in IVazza et al.i 
(t2011blh In this case, we differ from the analysis reported 
there by computing the azimuthal scatter from a larger 
number of angular sectors , N= 12, com pared to the cases of 
N=2, 4 and 8 explored in IVazza et al.l (l2011bli . In the sim- 
ulations, a number of dense clumps are present, which may 
bias the predicted scatter high. To overcome this problem, 
we computed the scatter of the simulated clusters both for 
the total gas distribution and by filtering out the 1% most 
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Fig. 9. Comparison between the average observed az- 
imuthal scatter profile from Fig. [7] (black) and the scatter in 
the simulations for the ENZO runs (red), for the total scat- 
ter (solid line) and when filtering out the 1% most-luminous 
cells (dashed curve). The cyan (non-radiative) and magenta 
(CSF) curves represent the scatter in the ART simulations. 


X-ray luminous cells, as in iRoncarelli et ah 12006 ) . which 
removes a large fraction of the clumps. 

We also performed a similar analysis on the set of ART 
simulations, both for the non-radiative and CSF runs. In 
this case, we analyzed mock X-ray images using the same 
method as the observational data (see Sect. 13.41) . and ap- 
plied our point-source detection algorithm to remove the 
most prominent clumps. Because of the small number of 
objects considered, we ignored the GADGET simulations for 
this analysis. For a comparison between GADGET and ENZO 
scatter profiles, we refer the reader to lVazza et al. 1 (12011b ). 

In Fig. [9] we show the measured scatter profile from 
Fig. El together with the scatter profiles of X-ray brightness 
from ENZO (red) and ART simulations (non-radiative: cyan; 
CSF: magenta). The dashed red profile shows the ENZO 
profile after filtering the 1% most-luminous pixels in the 
projected images. Interestingly, we note that non-radiative 
runs (red and cyan) overestimate the observed azimuthal 
scatter, while CSF simulations underestimate it. In the lat- 
ter case, radiative cooling is lowering the entropy of the 
gas, which makes it sink into the cluster’s potential well. 
This effect produces more spherical X-ray morphologies, 
thus lowering the azimuthal scatter. Conversely, in NR runs 
the effects of dynamics are more important, which create 
more substructures and increases the azimuthal scatter. 

Interestingly, the profile that best reproduces the data 
is the ENZO profile for which the 1% most-luminous pixels 
were filtered out (dashed red). This may indicate that some 
clumps are indeed present in the observations, but were 
detected as point sources and were masked for the analysis. 
We remark that even if in this case the azimuthal scatter 
from non-radiative simulation runs is in good agreement 
with the ROSAT data, the absolute profiles of density are 
too steep compared to observations (see the left panel of 
Fig. [HJ • However, our definition of the azimuthal scatter 
(Eq. [2]) is normalized to the absolute value of the profile at 


each radii, which makes it a rather robust proxy of cluster 
asymmetries over large ~ Mpc scales. 


6. Discussion 


6.1. Observational results 

In agreement w ith earlier works using R OSAT (V99, N05) 
and Chandra (lEttori fc Balest.ral 120091), but at variance 
with some recent resul ts from Suzaku ( Bautz et al.1 120091 : 
ISimiones cu et al. 2011 : George et al.1 120091) and XMM- 
Newton (lUrban et al.1 l201ll) . our analysis reveals that in 
average the slope of the density profiles steepens beyond 
’'soo (see Table H. This result indicates that the latter re- 
sults may have been performed along preferential directions 
connected with the large-scale structure (e.g., in the di- 
rection of filaments). Indeed, the narrow FOV of Suzaku 
allowed only a sparse coverage of the outskirts of nearby 
clusters, and hence these measurements might be the res ult 
of azim uthal variations. In the case of A1795. iBautz et ahl 
(I2009P detected a significant sign al only in the North di- 
rection, while the Perseus result (ISimionescu et al.l 120111) 
was obtained along two narrow arms, covering less than 
10% of the cluster’s extent at r 2 oo- Moreover, using sev- 
eral offset ROSAT /PSPC pointings of the Perseus cluster, 
lEttori et al.l (ll998l) observed clear azimuthal variations in 
density and gas fraction. Therefore, it is likely that the 
aforementioned measurements are not representative of the 
cluster as a whole. This picture is supported by our analysis 
of azimuthal variations in cluster outskirts, which suggests 
that even CC clusters exhibit significant departures from 
spherical symmetry around r 2 oo- Consequently, a full az- 
imuthal coverage is required to study the global behavior 
of cluster outer regions. 

An important and previously unknown result of this 
work is the systematic difference between CC and NCC 
cluster populations observed beyond ~ 0.3r2oo (see Fig. [4j) . 
As explained in Sect. [2 this effect seems to be an intrin- 
sic difference between the two classes, since it is does not 
correspond to a biased distribution of our sample in tem- 
perature or redshift. Our scaled gas mass profiles provide 
a natural explanation for this result (see Fig. [6]). Indeed, 
when the appropriate scaling is applied, the steeper den- 
sity profiles of CCs in the outskirts compensate exactly for 
the excess density in the central regions, such that clus- 
ters with the same virial mass have the same gas mass en- 
closed within r 2 oo, albeit distributed in a different way for 
relaxed and disturbed objects. This resu lt was expected in 
the old cooling- flow scenario (iFabianlll994l) , in which radia- 
tive cooling causes the gas to flow inwards and accumulate 
in the central regions. While in the central regions AGN 
feedback prevents the gas fr om cooling below a certain level 
(e.g.. (McNamara fe Nulsenl 12007P . the entropy injected by 
the central AGN is not sufficient to balance the flow in 
the outer regions of clusters, which explains the steep den- 
sity profiles seen in Fig. [I] Conversely, merging events are 
capable of injecting a very large amount of energy in the 
ICM, which results in an efficient redistribution of the gas 
between the core and the outer regions, and creates the 
flatter density profiles measured for NCC clusters. 

We also determined the typical scatter in surface- 
brightness as a function of radius (see Fig. and split 
the data into the CC and NCC classes. In the central re- 
gions, we observe a systematic difference between CC and 
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NCC clusters, NCC clusters showing a higher level of scat- 
ter than CC. This result is easily explained by the larger 
numb er of su bstru ctures gen erally observed in NCC clusters 
le.g.. ISanderson et al.ll2009t i. For CC clusters, we measure 
a scatter of 20 — 30% below 0.5r 2 oo, which corresponds to 
small variations (~ 10%) in gas density. This indicates that 
the azimuthal scatter in the inner regions (r < 0.5r2oo) 
can be used as an e stimator of the X-ra y state of clus- 
ters, as suggested bv IVazza et al.l (|2011bD . Conversely, in 
cluster outskirts the scatter of CC profiles increases, and 
there is no observed difference between the two classes. 
Interestingly, we note that for CC clusters the turnover in 
Fig. [3 occurs around rsoo, which coincides with the radius 
beyond which large scale infall motions and filame ntary 
accre tions are generally non- negligible (e.g., lEvrard et all 
ll99fii h Inside 7'5oo, the gas is virialized in the cluster’s po- 
tential well, and shows only little deviations from spherical 
symmetry. Beyond rsoo accretion processes are important, 
and the gas is located mostly along preferential directions 
(i.e. , filaments). As a result, the distribution of the gas be- 
comes strongly anisotropic, even for clusters which exhibit 
a relaxed morphology in their inner regions. 

6.2. Comparison with simulations 

Comparing our density profiles with numerical simulations, 
we find that all non-radiative simulations predict very steep 
profiles already starting from ~ 0 . 2 r 2 oo> with values of the 
/3 parameter larger than 0.85 in the 0.4 — 0.65r 2 oo range 
(see the left panel of Fig. [8] and Table (2). This indicates 
that the inclusion of non-gravitational effects is needed to 
reproduce the observed slope, even well outside of cluster 
cores. The runs including additional physics are in better 
qualitative agreement with the observations (see the right 
panel of Fig. 0, although because of overcooling their gas 
fraction is too low (~ 10% compared to ~ 15%). However, 
it seems unlikely that star formation and galactic winds (as 
in the CSF runs explored here) are the only necessary feed- 
back mechanisms needed to reproduce observed clusters. 
Indeed, simple feedback models still face severe problems 
in matching the properties of the stellar components inside 
galaxy clus ters, as w ell as the propert ies of galaxies within 
them (e.g.- lBorgani fo Kravtsovl[2009l for a recent review). 

As illustrated in Table [2 gas clumping may also play a 
role in reconciling simulations with observations. Indeed, if 
an important fraction of the gas in cluster outskirts is in 
t he for m of de nse g as clumps, as suggested in simulations 
(|Nagai fc haul l2CTTlh . the emissivity of the gas would be 
significantly increased, thus leading to an overestimation of 
the gas density when the assumption of constant density in 
each shell is made. Our results show that the treatment of 
gas clumping slightly improves the agreement between data 
and simulations (see the right panel of Fig. 0. In addition, 
gas clumping also provides an alternative interpretation for 
our observed difference between the CC and NCC popula- 
tions beyond 0.3r2oo- Indeed, simulations predict a larger 
clumping factor in unrelaxed clusters compared to relaxed 
systems for the same average density, which would result 
in a higher observed density in the former. At the moment, 
it is not clear whether this difference is caused by gas re- 
distribution or clumping, or if both of these effects play a 
role to some extent. 

On the other hand, we find that numerical simulations 
can reproduce qualitatively the observed azimuthal scat- 


ter in the galaxy cluster gas density profiles (see Fig. 0. 
Interestingly, we find that the observed azimuthal scat- 
ter is reproduced accurately when the 1% most-luminous 
clumps are filtered out, whereas the non-radiative simula- 
tions with no filtering overestimate the observed level of 
azimuthal scatter at all radii. Two possible interpretations 
can be put forward to interpret this result. Observationally, 
it is possible that the dense clumps were detected as point 
sources and were filtered out of our observations. If this 
is the case, long exposures with high-resolution X-ray tele- 
scopes ( Chandra or X MM-Newton) should allow us to char- 
acterize the point sources and discriminate between dense 
clumps and background AGN, possibly unveiling the popu- 
lation of accreting clumps in cluster outskirts. Conversely, 
if such observations do not confirm the existence of the 
clumps, it would imply that non-radiative simulations sig- 
nificantly over-estimate the amount of clumping in cluster 
outskirts, which would weaken the case for the interpreta- 
tion recently put forward to explain the fl attening of the 
entropy profiles observe d in a few cases (|Simionescu et al.l 
120111 : TlJrban et al]l2011li . 

As shown in Fig. 0 radiative cooling may also help 
to reconcile the non-radiative simulations with the data. 
Indeed, radiative cooling lowers the entropy of the gas and 
makes it sink into the potential well, which produces clus- 
ters with more spherical morphologies (iLau et al.ll2011f) and 
thus reduces the azimuthal scatter. Since we know that this 
effect is overestimated in our CSF simulations, radiative 
cooling likely reduces the azimuthal scatter with respect to 
non-radiative simulations, although not as much as what is 
predicted here. 

Alternatively, AGN feedback may be an important in- 
gredient which is rarel y taken int o acc ount in numerical 
simulations. Recently. iPratt et alJ (l2010f ) observed an anti- 
correlation between entropy and gas fraction, such that 
multiplying cluster entropy profiles by the local gas frac- 
tion allows to recover the entropy profiles predicted from 
adiabatic compression. I.e., the excess entropy observed in 
cluster cores is balanced by a lower gas fraction, and the to- 
tal entropy follows the p redictions of gravitational collapse. 
iMathews fc Guol (2011) interpreted this result in terms of 
the total feedback energy injected in the ICM through var- 
ious giant AGN outbursts, which they estimated to be as 
large as 10 63 ergs. In this scenario, feedback mechanisms are 
preventing the gas from collapsing into the potential well, 
causing a deficit of baryons in the inner regions of clusters, 
and thus flattening the observed density profiles. Moreover, 
it is well known that this mechanism also takes place at 
group and galaxy scale, leading to shallower density profiles 
in the accreting clumps. As a result, the gas distribution 
in cluster outskirts would be more homogeneous than pre- 
dicted in non-radiative simulations, in agreement with our 
observed azimuthal scatter profiles. Therefore, although its 
i mplementation int o numerical simulations is challenging 
(ISiiacki et al.ll2008lh AGN feedback could be an important 
effect to reconcile simulations with observations. A more 
complex picture of the ICM, possibly including also the de- 
tailed treatment of magnetic fields, cosmic rays, thermal 
conductions (and of the instabilities arising from these in- 
gredients), would still represent a challenge for present day 
cosmological simulations. 
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7. Conclusion 

In this paper, we presented our analysis of a sample of 
local (z = 0.04 — 0.2) clusters with i?OS^4T/PSPC, fo- 
cusing on the properties of the gas in cluster outskirts. 
We then compa red our observational results with numer- 
ical simulations (fRoncarelli et, al.ll2006l : iNagai fc Laull‘20 1 It 
IVazza et al.l 1201 lbh . Our main results can be summarized 
as follows. 


We observe a general trend of steepening in the ra- 
dial profiles of emission-measure and gas density be- 
yond ~ rsn o, in goo d agree m ent with e arlier works 
from IV ikhlinin et al. _J 1999H , iNeumani] (l2005lj and 
lEttori fc Balestraf i 2009 ) . As a result, the shallow den- 
sity profiles observed in several clusters by Suzaku 
dBautz et al.l 120091 : iSimionescu et al.l 1201 lh are proba- 
bly induced by observations in preferential directions 
(e.g., filaments) and do not reflect the typical behavior 
of cluster outer regions. 

We note for the first time a difference between the den- 
sity profiles of CC and NCC systems beyond ~ O.3r 20 o, 
which cannot be easily explained by any selection effect. 
We interpret this result by a different distribution of the 
gas in the two populations: the well-known density ex- 
cess in the core of CC clusters is balanced by a slightly 
steeper profile in the outskirts, which leads to the same 
gas mass enclosed within r 2 oo hr the two populations 
(see Fig. El) • Alternatively, this result could be caused 
by a larger clumping factor in disturbed objects, leading 
to an overestimate of the gas density of NCC clusters 
in the external regions. 

We also observe a clear difference in the azimuthal scat- 
ter between the two populations in the central regions, 
which is easily explained by the more disturbed mor- 
phology of NCC clusters. Conversely, beyond ~ r 5 oo 
both populations show a similar level of asymmetry (60- 
80%), which suggests that a signification fraction of the 
gas is in the form of accreting material from the large- 
scale structures. 

Comparing our ROSAT density profile with numeri- 
cal simulations, we find that all non-radiative numerical 
simulations fail to reproduce the observed shape of the 
density profile, predicting density profiles which are sig- 
nificantly too steep compared to the data (see Table [2] 
and Fig. [5). This implies that non-gravitational effects 
are important well outside the core region. The runs in- 
cluding additional physics (cooling, star formation, SN 
feedback) predict flatter profiles, although still too steep 
compared to the observations. Besides, it is well known 
that these simul ations over-pred ict the stellar fraction 
in clusters (|Borgani fe Kravtsov! l2009h . A slightly bet- 
ter agreement is found when a treatment of the observa- 
tional effects of gas clumping is adopted (INagai fc Lad 

IM1- 


— Non-radiative simulations are able to predict with good 
accuracy the observed azimuthal scatter profile, but 
only when the 1% most-luminous cells are filtered out 
(see Fig. [9J) . This result implies that either the clumps 
are quite bright and were masked as point sources in our 
analysis pipeline, in which case offset XMM-Newton and 
Chandra observations will be able to characterize them 
spatially and spectrally, or the non-radiative simulations 
significantly overestimate the effects of clumping on the 
observable X-ray properties. 


— As an alternative explanation, we suggest that AGN 
feedback might be important even at large radii, 
and could help to reconcile obse rvations an d sim- 
ulations. Indeed, rec ent works dPratt, et al.1 120101 : 
iMathews fe Guol 1201 111 indicate that feedback mecha- 
nisms may be responsible for the well-known deficit of 
baryons in cluster cores, thus leading to flatter gas dis- 
tributions out to large radii. Moreover, the existence of 
such mechanisms at group and galaxy scale could also 
dilute the accreting material at large radii, leading to a 
smaller azimuthal scatter. 
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Table 3. Master table of the cluster sample. Column descr iption: 1. Cluster name; 2. Effective exposure of the PSPC observation; 3. Redshift (from NED ): 4. 
Hydrogen column density, Nh, along the line of sight dKalberla et al.l [20051) : 5. Mean temperature in the 200-500 kpc radial range; 6. r^o fr om lArnaud et ahl 
(2005) scaling relations, in physical units; 7. Same as 6, in apparen t un its; 8. Central density n n (t his work); 9. Central entropy K n , from lCa,vagnolo et al.l (120091) : 
10. Reference for the temperature profile (l= ISnowden et all ( 2008| ): 2= ICavagnolo et all ((2009); 3= lDe Grandi fo Molendil (120021) ). 


Cluster 

Exposure [ks] 

z 

N h [10 22 cm" 2 ] 

fcTh 00—500 [keV] 

C 2 oo [kpc] 

r 200 [arcmin] 

no [10 3 cm 3 ] 

Kq [keV cm 2 ] 

Reference 

A85 

10.065 

0.05506 

0.028 

6.3 ±0.1 

1873 

29.17 

18.9 ±0.25 

12.5 

1 

A119 

14.758 

0.0442 

0.037 

5.0 ±0.1 

1673 

32.04 

2.1 ±0.34 

233.9 

2 

A133 

19.429 

0.0566 

0.0164 

4.0 ±0.09 

1494 

22.68 

14.0 ±0.18 

17.3 

1 

A401 

7.519 

0.07366 

0.0995 

7.9 ±0.15 

2077 

24.72 

5.3 ±0.66 

166.9 

2 

A478 

23.019 

0.0881 

0.0131 

6.56 ±0.08 

1883 

19.05 

18.8 ±0.19 

7.8 

1 

A644 

10.310 

0.0704 

0.0750 

7.7 ±0.1 

2054 

25.48 

9.4 ±0.29 

132.4 

2 

A665 

37.066 

0.1819 

0.0431 

8.0 ±0.2 

1987 

10.82 

5.6 ±0.18 

134.6 

1 

A1068 

10.822 

0.1375 

0.0173 

4.9 ±0.17 

1587 

10.89 

15.0 ±0.24 

9.1 

1 

A1651 

7.630 

0.084945 

0.0156 

6.7 ±0.2 

1913 

20.00 

8.8 ±0.50 

89.5 

2 

A1689 

14.291 

0.1832 

0.0186 

9.2 ±0.2 

2126 

11.51 

13.8 ±0.22 

78.4 

1 

A1795 

35.494 

0.06248 

0.0121 

6.02 ±0.08 

1828 

25.31 

20.1 ±0.12 

19.0 

1 

A1991 

21.956 

0.0586 

0.0248 

2.4 ±0.1 

1064 

15.64 

16.1 ±0.22 

1.5 

1 

A2029 

13.089 

0.07728 

0.0323 

7.7 ±0.2 

2054 

23.40 

20.2 ±0.20 

10.5 

1 

A2142 

19.410 

0.0909 

0.0383 

9.0 ±0.3 

2209 

21.73 

10.3 ±0.17 

68.1 

3 

A2163 

7.267 

0.203 

0.109 

18.8 ±1.3 

3008 

15.01 

8.2 ±0.92 

438.0 

2 

A2204 

5.346 

0.1526 

0.0561 

8.3 ±0.2 

2057 

12.93 

33.3 ± 0.76 

9.7 

1 

A2218 

43.179 

0.1756 

0.0266 

6.7 ±0.3 

1825 

10.22 

4.6 ±0.10 

288.6 

1 

A2255 

13.676 

0.0806 

0.0250 

6.1 ±0.1 

1817 

19.9 

2.3 ±0.32 

529.1 

2 

A2256 

17.000 

0.0581 

0.0418 

6.2 ±0.1 

1865 

27.63 

3.0 ±0.47 

349.6 

1 

A2597 

7.426 

0.0852 

0.0246 

3.64 ±0.06 

1405 

14.65 

18.0 ±0.22 

10.6 

1 

A3112 

7.829 

0.07525 

0.0137 

4.8 ±0.1 

1613 

18.82 

18.3 ±0.26 

11.4 

1 

A3158 

3.123 

0.0597 

0.0138 

5.1 ±0.1 

1681 

24.27 

3.8 ±0.20 

166.0 

1 

A3266 

13.967 

0.0589 

0.0158 

9.2 ±0.3 

2260 

33.05 

5.3 ±0.49 

72.5 

3 

A3558 

28.751 

0.048 

0.0402 

5.06 ±0.05 

1687 

29.89 

7.2 ±0.23 

126.2 

1 

A3562 

20.518 

0.049 

0.0376 

4.8 ±0.3 

1635 

28.41 

5.7 ±0.26 

77.4 

3 

A3667 

12.462 

0.0556 

0.0452 

5.31 ±0.05 

1721 

26.56 

4.5 ±0.36 

160.4 

2 

A4059 

5.684 

0.0475 

0.0122 

4.07 ±0.08 

1513 

27.08 

4.7 ±0.33 

7.1 

1 

Hydra A 

18.541 

0.0539 

0.0468 

4.0 ±0.06 

1495 

23.75 

22.1 ±0.17 

13.3 

1 

MKW 3s 

9.781 

0.045 

0.0272 

3.52 ±0.06 

1409 

26.54 

13.5 ±0.22 

23.9 

1 

PKS 0745-191 

9.627 

0.1028 

0.405 

8.4 ±0.3 

2121 

18.70 

31.9 ±0.45 

12.4 

1 

Triangulum 

7.343 

0.051 

0.114 

8.9 ±0.2 

2229 

37.31 

5.9 ±0.79 

313.0 

1 
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Appendix A: Notes on individual objects 

- A85: 

A sub-cluster located ~ 10' South of the cluster center 
is currently merging with the main cluster. This sub- 
structure was masked for the analysis. 

- A401: 

The cluster is connected through a filament to its neigh- 
bor A399, located ~ 35' South-West of the center of 
A401. We extracted the surface-brightness profile in a 
sector of position angle 340-250° to avoid any contami- 
nation of A399 to our measurement of the CXB. 

- A478: 

The combination of a favorable temperature/redshift 
and of a good-quality ROSAT observation allow us to 
reach the highest signal-to-noise ratio in the sample at 
r 2 oo for this strong CC cluster. As a result, the data from 
this cluster may contribute strongly when a weighted 
mean is performed. 

- A644: 

This NCC cluster exhibits an unusual decreasing az- 
imuthal scatter profile, showing a large (close to 100%) 
scatter in its central regions, but no significant scatter 
around r 2 oo- 

- A 2029: 

A probable filament is connecting A2029 to A2033, lo- 
cated ~ 35' North of the center of A2029. The surface- 
brightness profile was extracted in a sector with position 
angle 140-80° to measure the CXB level. 

- A2142: 

Several PSPC observations of this famous cold-front 
cluster exist. For this work, we used the longest avail- 
able observation, which was pointed 16' South of the 
center of A2142. This is the only case in the sample for 
which the observation was not pointed on the target. 

- A3 5 58 and A3 5 62: 

These two clusters are located in the Sliapley super- 
cluster and connected by a filament. Consequently, they 
show an unusually high azimuthal scatter in the out- 
skirts. The CXB level was estimated by excluding the 
direction of the filament. 

- A3 667: 

This very disturbed cluster shows the highest emission- 
measure and density in the sample beyond ~ 0.2r 2 oo, 
and hence it could bias our average profiles, in particu- 
lar when computing the difference between the CC and 
NCC classes. However, removing it from the sample did 
not lead to any significant difference, either quantitative 
or qualitative. 

- A4059: 

This is the most azimuthally-symmetric cluster in the 
sample. The azimuthal scatter for this cluster is consis- 
tent with 0 at all radii. 

- Hydra A: 

A tail of emission (filament?) extends out to ~ 20' 
South-East of the cluster core. This leads to a very high 
azimuthal scatter (> 100%) around r 20 o- 


Appendix B: Mean emission-measure profiles 

In Table EH we give the mean self-similar scaled 
emission-measure profiles for the CC and NCC classes 
and the whole sample, as shown in Fig. Q] 
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Table B.l. Data of Fig. QJ mean self-similar scaled emission-measure profiles for the whole sample and for the CC and 
NCC classes, in units of cm ~ 6 Mpc 


Rin 

Rout 

Total 


CC 


NCC 


0 

0.02 

(1.80 ±0.01) • 

10 _s 

(9.48 ± 0.05) ■ 

10 _s 

(1.13 ±0.01) 

• 10 _t> 

0.02 

0.04 

(1.26 ±0.01)) 

10~ 5 

(4.83 ± 0.02) • 

10“ 5 

(8.32 ± 0.06) 

• 10" 6 

0.04 

0.06 

(9.63 ± 0.04) ■ 

10" 6 

(2.28 ±0.01) • 

10“ 5 

(6.90 ± 0.04) 

• IO" 6 

0.06 

0.08 

(7.39 ±0.03) ■ 

10" 6 

(1.23 ±0.01) • 

10“ 5 

(5.70 ±0.03) 

• 10~ 6 

0.08 

0.1 

(5.45 ± 0.02) ■ 

10" 6 

(7.72 ± 0.04) ■ 

10" 6 

(4.49 ± 0.02) 

• 10" 6 

0.1 

0.12 

(4.12 ±0.02) • 

10" 6 

(5.27 ±0.03) • 

10" 6 

(3.52 ± 0.02) 

• IO" 6 

0.12 

0.14 

(3.20 ± 1.36) ■ 

10" 6 

(3.63 ±0.02) • 

10" 6 

(2.91 ±0.02) 

• IO" 6 

0.14 

0.16 

(2.47 ±0.01) ■ 

10" 6 

(2.60 ±0.02) • 

10" 6 

(2.37 ±0.01) 

• 10' 6 

0.16 

0.18 

(1.91 ±0.01) ■ 

10" 6 

(1.95 ±0.01) • 

10" 6 

(1.88 ±0.01) 

• 10~ 6 

0.18 

0.2 

(1.51 ±0.01) • 

10" 6 

(1.48 ±0.01) • 

10" 6 

(1.54 ±0.01) 

• 10" 6 

0.2 

0.22 

(1.23 ±0.01) • 

10" 6 

(1.19 ±0.01) • 

10" 6 

(1.26 ±0.01) 

• IO -6 

0.22 

0.24 

(1.02 ±0.01) • 

l 0 -6 

(9.47 ±0.09) ■ 

IO" 7 

(1.07 ±0.01) 

• IO" 6 

0.24 

0.26 

(8.40 ± 0.05) ■ 

10“ 7 

(7.61 ±0.08) ■ 

IO" 7 

(8.95 ± 0.07) 

• 10“ 7 

0.26 

0.29 

(6.91 ±0.05) ■ 

10" 7 

(6.09 ± 0.07) • 

IO" 7 

(7.59 ± 0.06) 

• 10“ 7 

0.29 

0.31 

(5.32 ± 0.04) • 

10“ 7 

(4.73 ± 0.06) • 

10“ 7 

(5.77 ±0.05) 

• 10“ 7 

0.31 

0.34 

(4.30 ± 0.04) • 

10“ 7 

(3.74 ± 0.06) ■ 

10“ 7 

(4.70 ± 0.05) 

• 10~ 7 

0.34 

0.38 

(3.20 ± 0.03) ■ 

10“ 7 

(2.77 ±0.04) ■ 

10“ 7 

(3.60 ± 0.04) 

• IO" 7 

0.38 

0.41 

(2.49 ± 0.02) ■ 

10" 7 

(2.10 ±0.04) ■ 

10“ 7 

(2.76 ±0.03) 

• 10~ 7 

0.41 

0.45 

(1.86 ±0.02) ■ 

10“ 7 

(1.57 ±0.03) • 

10“ 7 

(2.11 ±0.03) 

• 10' 7 

0.45 

0.50 

(1.48 ±0.02) ■ 

io- 7 

(1.27 ±0.03) • 

io- 7 

(1.63 ±0.02) 

■ io- 7 

0.50 

0.55 

(1.07 ±0.02) ■ 

10“ 7 

(9.05 ± 0.24) • 

10 -8 

(1.18 ±0.02) 

• 10“ 7 

0.55 

0.60 

(7.99 ±0.14) ■ 

10 -8 

(6.82 ±0.22) • 

10 -8 

(8.87 ±0.19) 

• 10" 8 

0.60 

0.66 

(5.73 ±0.12) ■ 

10" 8 

(4.97 ±0.18) • 

10" 8 

(6.30 ±0.16) 

• 10 -8 

0.66 

0.72 

(4.28 ±0.11) ■ 

10" 8 

(3.78 ± 0.17) • 

10" 8 

(4.62 ±0.14) 

• 10" 8 

0.72 

0.79 

(3.06 ±0.11) ■ 

10" 8 

(2.75 ±0.18) • 

10" 8 

(3.21 ±0.13) 

• 10 -8 

0.79 

0.87 

(2.23 ±0.10) ■ 

10 -8 

(1.77 ±0.16) • 

10 -8 

(2.51 ±0.13) 

• 10" 8 

0.87 

0.95 

(1.35 ±0.09) • 

10 -8 

(8.57 ± 1.49) • 

10“ 9 

(1.63 ±0.11) 

• 10' 8 

0.95 

1.05 

(7.77 ± 0.85) • 

10" 8 

(5.85 ± 1.40) • 

IO” 9 

(8.88 ± 1.07) 

• 10~ 9 

1.05 

1.15 

(5.32 ±0.80) ■ 

10~ 8 

(4.19 ± 1.35) • 

10" 9 

(5.92 ±0.99) 

• 10 -9 

1.15 

1.26 

(4.74 ±0.81) • 

10" 8 

(3.75 ± 1.40) • 

10" 9 

(5.24 ±0.97) 

• 10 -9 


